Interlaced Dynamical Decoupling and Coherent Operation of a Singlet-Triplet Qubit 
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We experimentally demonstrate coherence recovery of singlet-triplet superpositions by interlacing 
qubit rotations between Carr-Purcell (CP) echo sequences. We then compare performance of Hahn, 
CP, concatenated dynamical decoupling (CDD) and Uhrig dynamical decoupling (UDD) for singlet 
recovery. In the present case, where gate noise and drift combined with spatially varying hyperfine 
coupling contribute significantly to dephasing, and pulses have limited bandwidth, CP and CDD 
yield comparable results, with T^ ~ 80 [is. 

PACS numbers: 



The singlet (S = (\U) ~ lit)) /V$) and m = triplet 
(T = (|tl) + lit)) /V%) spin states of two electrons in 
a double quantum dot form a versatile qubit that is in- 
herently protected from collective dephasing pQ and al- 
lows sub-nanosecond rotation around one axis via elec- 
trical control of the exchange interaction [2]. Rotations 
around a second axis, needed for universal control, can 
be induced by a Zeeman field difference between the two 
quantum dots, created by a proximal micromagnet [3] or 
controlled Overhauser fields [4|. Fast measurement of a 
qubit state has been demonstrated for two-electron spin 
states [5HZ] , and also for single electron spin states [8j [9] . 

Temporal fluctuations of Overhauser fields evolve 
slowly compared to gate operation times [10 , allowing 
a simple one-pulse Hahn echo jTT] to significantly extend 
qubit coherence [2] . By repeating 7r-pulses on shorter in- 
tervals, Carr-Purcell (CP) sequences [I2j [13] can extend 
qubit coherence to impressively long times [14]. Theo- 
retical work and experiments in other systems suggests 
that, depending on the environmental noise spectrum, 
echo sequences more complex than CP may further in- 
crease coherence times [15-20 ]. To date, experimental 
demonstration of decoupling schemes for spin qubits have 
been limited to recovery of an initially prepared singlet 
state. For applications in quantum information process- 
ing, however, decoupling schemes must preserve an arbi- 
trary qubit state [21]. 

In this Letter, we demonstrate echo recovery of singlet- 
triplet superposition amplitude by embedding qubit ro- 
tations about two axes of the Bloch sphere between CP 
sequences. For Overhauser-driven rotations about the x 
axis of the Bloch sphere, and exchange- driven rotations 
about the z axis, ensembles of single-shot singlet mea- 
surements show coherent oscillations for total sequence 
times (including echoes) of 60 /is. We then compare sin- 
glet recovery using Hahn, CP, concatenated dynamical 
decoupling (CDD) [15J[T6], and Uhrig dynamical decou- 
pling (UDD) [17] schemes. We find that for the present 
setup, a 16-pulse CP sequence and 21-pulse (fifth order) 
CDD sequence yield comparable performance, both with 
T2 ~ 80 /is, while a 22-pulse (22nd order) UDD performs 



less well, possibly due to the limited bandwidth of the 
pulses. 

The double quantum dot and sensor are defined by 
Ti/Au depletion gates on a GaAs/Al .3Gao.7As het- 
erostructure with a two-dimensional electron gas (den- 
sity 2 x 10 15 m -2 , mobility 20 m 2 /Vs) 100 nm below the 
surface. An in-plane magnetic field, B = 750 mT, was 




FIG. 1: (Color online) (a) Micrograph of lithographically 
identical device with contacts, dot locations, and GaAs crys- 
tal axes indicated. Gate voltages, VL and Vr, set the energy 
of left and right dot. An rf-coupled sensor dot measures the 
double dot charge state, (b) dc conductance change, Ag, 
with double dot charge state. Markers indicate gate voltages 
during experimental pulse sequences. Detuning, e, is con- 
trolled by Vl and Vr along the diagonal through points S 
and M. (c) Energy level diagram as function of detuning e, 
indicating locations of ep where the (0,2) singlet is prepared, 
es of separated electrons for dephasing and x-rotations, ce 
of exchange pulse and 7r-pulses, 6r of mapping ramps, and 
cm of the measurement point. Exchange energy, Je used for 
7T pulses is indicated, (d) Bloch sphere of the singlet-triplet 
qubit, with mechanisms of rotation indicated, exchange en- 
ergy, J, and Zeeman field difference along the applied field 
direction, AB Z , drive rotations around the qubit z-axis and 
x-axis respectively. 



applied perpendicular to the dot connection axis, as in- 
dicated in Fig. [TJa). Measurements were carried out in 
a dilution refrigerator (electron temperature ~ 150 mK) 
configured for high- bandwidth gating and rf reflectome- 
try. The state of the double quantum dot was controlled 
by pulsed gate voltages Vl, Vr [Figs.fl[a,b)] using a Tek- 
tronix AWG 5014. Except for point P, where the (0,2) 
singlet is prepared, gate configurations fall on a line be- 
tween (0,2) and (1,1) charge states, parameterized by the 
detuning, e [Fig. [lib)]. To reduce effects of voltage drift, 
average gate voltages were set to the separation point S 
using a compensation pulse between measurement and 
preparation [14] . The charge state of the double dot was 
detected using rf reflectometry [22] of a proximal sensor 
quantum dot, integrated over ~ 600 ns to yield a single- 
shot measurement, as described previously [6j[7]- 

All pulse sequences start by allowing two electrons at 
point P to relax to the (0,2) singlet [Fig. 1(c)]. Fast 
separation of the electrons to point S initializes the sys- 
tem into the (1,1) singlet, the +z direction on the qubit 
Bloch sphere [Fig. 1(d)]. Alternatively, fast separation 
from P to R followed by adiabatic separation from R 
to S initializes the system into the (1,1) ground state 
of the Overhauser fields, |t|), the +x direction on the 
Bloch sphere [2]. At point S, exchange splitting, J$, 
is negligible but Overhauser gradients are not, and the 
qubit precesses rapidly around its x axis at frequency 
/ s = g*/j jB AB z /h, where g* ~ -0.4 is the GaAs g- 
factor and AB Z is difference in Zeeman fields (in this 
case, Overhauser fields) along the direction of applied 
field. Controlled rotation about the z-axis of the Bloch 
sphere is realized by pulsing from S to E, where a large 
exchange energy, Je, causes qubit rotation at frequency 
/e = y/Jn + (g*fiBAB z ) 2 /h. Frequency f E includes a 
small quadrature contribution from AB Z that causes in- 
homogeneous dephasing and alters the rotation axis dur- 
ing the exchange pulse [Fig. 1(d)]. After evolution at 
S and E, the qubit is measured by pulsing to the mea- 
surement point, M. The singlet state can recombine in 
(0,2) while the triplet state remains in (1,1). The re- 
sulting charge state difference is detected by the sensor 
quantum dot as V T { [6j [7] . 

Dephasing of the S-Tq qubit can arise from thermally 
driven evolution of the Overhauser field difference be- 
tween dots or from electron motion in a spatially varying 
Overhauser field caused by gate noise or drift. Here, we 
do not distinguish between decoherence due to quantum 
entanglement with the environment and dephasing due 
to a noisy classical environment. Because Overhauser 
fluctuations are concentrated below 1 Hz [10 , apparent 
dephasing due to precession of the qubit can be readily 
recovered using dynamical decoupling schemes applied on 
faster time scales [2j[l4]. The simplest such scheme, Hahn 
echo (HE), uses a single 7r-pulse, realized in this case by 
pulsing to E for a time tt/^HJe) after a (dephasing) inter- 
val td/2 at S, followed by a second (rephasing) interval 
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FIG. 2: (Color online) (a) Pulse sequence for x rotations, cy- 
cling through preparation (P), separation (S), exchange-echo 
(E), and measurement (M) points in gate- voltage space. Evo- 
lution at S for time rs produces x rotations due to Overhauser 
gradient, (b) Singlet return probability, Ps, as function of rs 
for Ntt — 4 (green diamonds) and N n = (black circles), 
with td = 4 /xs, along with cosine fit (curve) to data, post- 
selected to have a specific period (see text). Without echo 
pulses (Nn — 0), no oscillations are seen. Error bar reflects 
statistics of 100 single-shot measurements per point, (c) Iden- 
tical to (b) with td = 16 /xs, showing reduced visibility, (d) 
Pulse sequence for z rotations: after an initial CP sequence, 
singlet is mapped to |ti) by pulsing to R and ramping to S. 
Pulsing to E for te induces rotation around z. Pulsing to S 
then ramping to R then pulsing back to S maps |ti) back 
into the singlet. After a second CP series, pulsing to M yields 
single-shot measurement, (e) Singlet probability, Ps, as func- 
tion of rs for Nn = 8 (blue diamonds) and N n = (black 
circles), with td = 16 /xs, based on an average of 600 x 100 
single-shot measurements, along with cosine fit with gaussian 
envelope, see text. Error bars from statistics are smaller than 
markers, (f) Identical to (e) with td = 64 /xs. 



at S for td/2. Multiple 7r-pulses applied within a given 
interval can extend coherence by shortening the inter- 
val between dephasing and rephasing, during which the 
environment may have evolved [12] . The simplest multi- 
pulse scheme, the CP sequence, inserts several of equally 
spaced 7r-pulses surrounded by dephasing and rephasing 
intervals [T2lfT3] . 

As a first approach to dynamical decoupling of arbi- 
trary qubit states, we investigate x and z rotations in- 
terlaced between two CP sequences. For this purpose, 
the sequence preceding the rotations is not needed, but 



including it serves to demonstrate general interlacing of 
qubit operations with dynamical decoupling. 

To interlace x rotations between CP sequences [see 
Fig. [2|a)], the (1,1) singlet was first initalized, followed 
by a N n = 4 CP sequence during an interval td/2. Paus- 
ing at S for a time r$ induced rotation around x axis 
due to the slowly- varying Overhauser field gradient. This 
was followed by a second N n = 4 CP sequence, and fi- 
nally a single-shot measurement at point M. For each 
value of rs, this sequence was repeated 100 times, yield- 
ing a probability, Ps, of a singlet outcome over these 
100 single-shot events. Probabilities for 40 values of rs, 
ranging from to 180 ns, were measured over a total du- 
ration of 100 ms. Oscillations of Ps(ts) are well fit by 
Ps(r s ) = 0.5 (1 + Vcos (r s ^/i B A5 z //i)), where V is the 
measurement visibility [Figs. 2(b,c)], the expected form 
for qubit precession in a fixed Overhauser field difference 
AB Z [6]. Averaging x rotations over the full range of 
Overhauser fields would wash out oscillations. Instead, 
Ps(t~s) data in Figs. 2(b,c) was selected to cover a nar- 
row range of Overhauser fields around AB Z ~ 3 mT. The 
value of AB Z evolves slowly among data sets, observable 
as a fluctuating oscillation period of Ps(ts). 

From the ratio of visibilities for long and short (but 
otherwise identical) sequences, V(jd = 16/is) /V(td = 
4/is) = 0.88, we estimate a dephasing time T 2 ~ 40 /is for 
the complete sequence with N n = 4 before and after the 
rotation, assuming a gaussian decay envelope (discussed 
below) and using the measured normalizations Vo ~ 0.70 
with r D = 0, Nn = 0, and V ~ 0.63 with r D = 0, N n = 4. 
The reduced visibility of S-Tq precession without the CP 
sequences can be attributed to limited read-out fidelity, 
finite exchange at es, and imperfect initialization in (1,1). 
A further ~ 10% reduction in visibility when including 
2N n = 8 pulses is dominated by tilting of the echo ro- 
tation axis away from z by Overhauser field gradients 
[Fig.Qd)]. 

Interlaced z rotations between CP sequences were 
demonstrated using adiabatic ramps that map S and To 
to |tl) and ||t) [Fig.^d)] 0. Initialization of the (1,1) 
singlet was followed by a N n = 8 CP sequence. The pre- 
served singlet was then mapped to |t|) by pulsing from S 
to R (without passing the S-T+ anticrossing) then ramp- 
ing slowly to S in a (loading) time tl ~ 400ns. Pulsing to 
E for a time te induced rotation around z at frequency 
/e- Next, pulsing to S, ramping slowly from S to R, then 
back to S mapped |tl)(|lt)) to S(T ). The resulting su- 
perposition amplitude was preserved over a time td/2 by 
a second N n = 8 CP sequence. Finally, the superposi- 
tion amplitude was measured by pulsing to M. 600 sin- 
glet probabilities, each based on 100 single-shot measure- 
ments, were averaged over fluctuating Overhauser fields 
for each te- Oscillations in Ps(te) do survive averaging 
over Overhauser fields, since for Je ^> AB Z the oscil- 
lation frequency /e depends only weakly on AB Z . How- 
ever, this averaging does contribute to the decay envelope 
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FIG. 3: (Color online) (a) Singlet recovery amplitude, Ps, as 
a function of total dephasing time, td, for Hahn, CP, CDD 
and UDD sequences (see text), with fits of Eq. (HI with a — 2, 
yielding T 2 HE ~ 6 /is, T 2 CDD ~ 80 /is, T 2 CP ~ 80 /is, T 2 UDD ~ 
30 /is ^26 . (b) Replotting data in (a) as -ln[2(P s - 1/2) /V] 
allows a fit giving the exponent a in Eq. (1). On a log-log 
plot, Eq. (1) and data appears linear with slope a [26]. Pulse 
sequences for (c) 5th order CDD and (d) 22nd order UDD. 



seen in Figs. 2(e,f). Note that mapping between z and 
x axes accurately preserves the angle of the superposi- 
tion from the corresponding axis (amplitude) , with some 
loss of azimuthal angle (phase) information due to envi- 
ronment dynamics on the time scale of the 400 ns ramp. 
State tomography [4 requires fast rotations compared to 
environment dynamics [TO] . 

Singlet return probability Ps'(te), averaged over fluc- 
tuating Overhauser fields, for td = 16 /is [Fig. 2(e)] and 
td = 64 /is [Fig. 2(f)] using N n = 8 CP sequences show 
damped oscillations a function of the exchange time te- 
Fits to P s = 111 + (F/2)cos(rEJEA)exp[-(r E /T 2 *) 2 ], 
the expected form for a qubit precessing in exchange 
splitting, Je, with inhomogenous dephasing due to both 
fluctuating Overhauser fields and electrical noise, yield 
Je ~ 0.2 fieV and T 2 * ~ 100 ns from the decay enve- 
lope [23 . Visibilities V - 0.09(0.06) for r D = 16(64) /is, 
along with V ~ 0.3 measured without CP sequences (not 
shown), yield a coherence time T2 ~ 80 /is. Besides fi- 
nite coherence time (T2), other contributions that reduce 
visibility include nonzero AB Z during exchange pulses, 
nonzero exchange at the separation point, Js, and gate 
noise during 7r-pulses. 

Next, we consider alternative dynamical decoupling 
schemes for recovery of a prepared singlet without in- 
terlaced rotations. A prepared (1,1) singlet at point S 
was given n^ 7r-pulses, spaced in time according to the 
particular decoupling scheme, over a total interval td. 
For each value of td, an Overhauser- averaged singlet re- 
turn probability Ps was found by averaging ~ 10 4 single- 
shot measurements, normalized by reflectometer output 



voltages of singlet and triplet outcomes via single-shot 
histograms |B]. For all decoupling schemes, Ps(tb) data 
are well described by the functional form 



P s (t d ) = 1/2 + (V/2)e 



-(td/T 2 )° 



(1) 



We either set a = 2 [Fig. 3(a)] or leave a as a fit pa- 
rameter [Fig. 3(b)]. Allowing a to vary can give insight 
into the dominant dephasing mechanism, depending on 
the decoupling scheme [24j [25] . For instance, for sim- 
ple Hahn echo, a white- noise environment is expected to 
give exponential decay, a — 1. In contrast, the enhanced 
low-frequency content of an Overhauser-field-dominated 
environment is expected to yield a ~ 4 for Hahn echo 
[24| [25] , as observed experimentally [14] . 

Experimental Ps(td), for several decoupling schemes 
along with fits of Eq. (fil) with fixed a = 2 yields val- 
ues for T2 [Fig. 3(a)]. For a single-pulse Hahn echo, 
T 2 HE - 6 /is, while T 2 CP - 80 /is for a CP sequence 
with 16 7r-pulses [26 . A comparison of CP sequences 
with different N n are shown in Fig. 3(b). Best-fit values 
for a are extracted by taking the logarithm of Eq. (1), 
as shown in Fig. p[c). For Hahn echo, the best fit value 
is a RE = 1.8, suggesting that electrical (gate) noise and 
drift, combined with spatially varying Overhauser fields, 
are likely the dominant source of dephasing, rather than 
time dependence of the Overhauser fields themselves. For 
the CP sequence, the best fit exponent is a CP ~ 1.5. 

Alternative decoupling sequences may outperform 
CP, depending on the decohering environment [T6HT8] . 
A favorable scheme for spin qubits coupled to nuclear 
environments is CDD whose nth order sequence is 
created from two sequences of (n — l)th order with 
an additional 7r-pulse between them for odd n [T6] : 
first-order CDD is Hahn echo; second-order CDD is CP 
with two 7r-pulses; fifth-order CDD, with 21 7r-pulses, is 
shown in Fig. [3lc). Another pulse scheme optimized for 
spin-bath environments is UDD, which has n^ 7r-pulses 
(indexed by j) at times Sr = td sin 2 [j7r/(2n + 2)]. 
The 22nd-order UDD (n^ = 22) is shown in Fig^d). 
Singlet probabilities Ps(td) for the 5th-order CDD 
and 22nd-order UDD sequence is shown in Fig. |3ja) 
along with Hahn and CP, and yield T 2 CDD ~ 80 /is and 
T 2 UDD - 30 /is [26 . We note that (i) all multi-pulse 
sequences significantly outperform Hahn echo; (ii) the 
fifth-order CDD sequence (n n = 21) has no better 
performance than the n n = 16 CP sequence, and 
(iii) the T^ achieved for UDD is considerably shorter 
than for comparable CDD and CP sequences. This 
presumably results from experimental artifacts such 
as pulse bandwidth limitations, but could also reflect 
the non-optimality of UDD for an environment with a 
1/f 2 high-frequency tail [10 J. For a dephasing power 
spectrum closer to the experiment the optimal spin echo 



pulse sequence has been predicted to be very similar to a 
simple CP sequence [20], while in Ref.[T8lit was predicted 
that UDD is optimal for hyp erfine- induced dephasing 
in GaAs quantum dots. Optimizing pulse sequences 
with contributions from gate noise and drift combine 
with spatially and temporally varying Overhauser fields 
remains an outstanding problem experimentally and 
theoretically. 
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